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Abstract 

Here we summarize the results from our direct numerical simulations (DNS) and experimental measurements on 
rotating Rayleigh-Benard (RB) convection. Our experiments and simulations are performed in a cylindrical samples 
with aspect ratio of 0.5 < T < 2.0. Here T = D/L with D and L are the diameter and height of the sample, respectively. 
When the rotation rate is increased, while a fixed temperature difference between the hot bottom and cold top plate 
is maintained, a sharp increase in the heat transfer is observed before the heat transfer drops drastically at stronger 
rotation rates. Here we focus on the question of how the heat transfer enhancement with respect to the non-rotating 
case depends on the Rayleigh number Ra, the Prandtl number Pr, and the rotation rate, indicated by the Rossby 
number Ro. Special attention will be given to influence of the aspect ratio on rotation rate that is required to get heat 
transport enhancement. In addition, we will discuss the relation between the heat transfer and the large scale flow 
structures that are formed in the different regimes of rotating RB convection and how the different regimes can be 
identified in experiments and simulations. 

Keywords: Rotating Rayleigh-Benard convection, direct numerical simulation, experiments, heat transfer, different 
turbulent states, flow structure 



Rayleigh-Benard (RB) convection, i.e. the flow of 
a fluid heated from below and cooled from above, is 
the classical system to study thermally driven turbu- 
lence in confined space (HE). Buoyancy-driven flows 
play a role in many natural phenomena and techno- 
logical applications. In many cases the fluid flow is 
also affected by rotation, for example, in geophysical 
flows, astrophysical flows, and flows in technology 0. 
On Earth, many large-scale fluid motions are driven by 
temperature-induced buoyancy, while the length scales 
of these phenomena are large enough to be influenced 
by the Earth's rotation. Key examples include the con- 
vection in the atmosphere (4) and oceans |5 ], including 
the global thermohaline circulation (51 . These natural 
phenomena are crucial for the Earth's climate. Rotating 
thermal convection also plays a significant role in the 
spontaneous reversals of the Earth's magnetic field Q. 
Rotating RB convection is the relevant model to study 
the fundamental influence of rotation on thermal con- 
vection in order to better understand the basic physics 
of these problems. 
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In this paper we discuss the recent progress that has 
been made in the field of rotating RB convection. First 
we discuss the dimensionless parameters that are used to 
describe the system. Subsequently, we give an overview 
of the parameter regimes in which the heat transport in 
rotating RB is measured in experiments and direct nu- 
merical simulations (DNSs). This will be followed by a 
description of the characteristics of the Nusselt number 
measurements and a description of the flow structures 
in the different regimes of rotating RB. Finally, we ad- 
dress how the different turbulent states are identified in 
experiments and simulations by flow visualization, de- 
tection of vortices, and from sidewall temperature mea- 
surements. 



1. Rotating RB convection 

When a classical RB sample is rotated around its cen- 
ter axis, it is called rotating RB convection. For not 
too large temperature gradients, this system can be de- 
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scribed with the Boussinesq approximation 
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for the velocity field u, the kinematic pressure field p, 
and the temperature field 6 relative to some reference 
temperature. In the Boussinesq approximation it is as- 
sumed that the material properties of the fluid such as 
the thermal expansion coefficient (3, the viscosity v, and 
the thermal diffusivity k do not depend on tempera- 
ture. Here g is the gravitational acceleration and O is 
the rotation rate of the system around the center axis, 
pointing against gravity, O = Qz. Note that p(7o) has 
been absorbed in the pressure term. In addition to the 
Oberbeck-Boussinesq equations we have the continuity 
equation (V-u = 0). A no-slip (u = 0) velocity boundary 
condition is assumed at all the walls. Moreover, the hot 
bottom and the cold top plate are at a constant tempera- 
ture, and no lateral heat flow is allowed at the sidewalls. 

Within the Oberbeck-Boussinesq approximation and 
for a given cell geometry, the dynamics of the system 
is determined by two dimensionless control parameters, 
namely, the Rayleigh number 



Ra = 



PgAL 3 
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where L is the height of the sample, A = 7^ - T t the 
difference between the imposed temperatures 7& and T t 
at the bottom and the top of the sample, respectively, 
and the Prandtl number 



Pr= V -. 

K 
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The rotation rate is indicated by the Rossby number 
Ro = V^A/L/(2i2), 



(5) 



which indicates the ratio between the buoyancy and 
Coriolis force. Note that the Ro number is an inverse 
rotation rate. Alternative parameters to indicate the ro- 
tation rate of the system are the Taylor number 
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comparing Coriolis and viscous forces, or the Ekman 
number 

Ek=-^ = ^=. (7) 

A convenient relationship between the different dimen- 
sionless rotation rates is Ro = ^Raf(PrTa). 



The cell geometry is described by its shape and an 
aspect ratio 

r = diu (8) 

where D is the cell diameter. The response of the sys- 
tem is given by the non-dimensional heat flux, i.e. the 
Nusselt number 
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where Q is the heat-current density and A the thermal 
conductivity of the fluid in the absence of convection, 
and a Reynolds number 



Re = 



UL 



(10) 



There are various reasonable possibilities to choose a 
velocity U, e.g., the components or the magnitude of 
the velocity field at different positions, local or averaged 
amplitudes, etc, and several choices have been made by 
different authors. 



2. Parameter regimes covered 

In figure [T] we present the explored Ra - Pr - Ro 
parameter space for rotating RB convection Here 
we emphasize that numerical simulations and experi- 
ments on rotating RB convection are complementary, 
because different aspects of the problem can be ad- 
dressed. Namely, in accurate experimental measure- 
ments of the heat transfer a completely isolated sys- 
tem is needed. Therefore, one cannot visualize the flow 
while the heat transfer is measured. On the positive side, 
in experiments one can obtain very high Ra numbers 
and long time averaging. In simulations, on the other 
hand, one can simultaneously measure the heat trans- 
fer while the complete flow field is available for anal- 
ysis. But the Ra number that can be obtained in sim- 
ulations is lower than in experiments, due to the com- 
putational power that is needed to fully resolve the flow. 
Here we should also mention that the highest Ra number 
reached in rotating RB experiments is almost Ra = 10 16 , 
whereas in DNSs of rotating RB convection the highest 
Ra number is Ra = 4.52 x 10 9 . However, the flexibil- 
ity of simulations allows to study more Pr numbers, i.e. 
covering a range of 0.2 < Pr < 100, whereas present ex- 
periments are almost exclusively for 3.05 < Pr < 6.4, 
i.e. the Pr number regime accessible with water. In ad- 
dition, we note that the 1 /Ro number regime that can 



figure 1 of Stevens et al. QD also the Ra-Ro-T for Pr = 4.38 
can be found. 
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be covered in experiments can, depending on Ra and 
Pr, be somewhat limited. Very low 1/Ro values, corre- 
sponding to very weak rotation, are difficult due to the 
accuracy limitations at very small rotation rates. The 
highest l/Ro that can be obtained in a setup is either 
determined by the requirement that the Froude number 
Fr = Q?(L/2)/g (9), which indicates the importance of 
centrifugal effects, is not too high (usually Fr < 0.05 is 
considered to be a safe threshold ifTOli ) or by the maxi- 
mum rotation rate that can be achieved. 



3. Nusselt number measurements 

Early linear stability analysis, see e.g. Chandrasekhar 
|271, revealed that rotation has a stabilizing effect due to 
which the onset of heat transfer is delayed. This can be 
understood from the thermal wind balance, which im- 
plies that convective heat transport parallel to the ro- 
tation axis is suppressed. Experimental and numeri- 
cal investigations concerning the onset of convective 
heat transfer and the pattern formation in cylindrical 
cells just above the onset under the influence of rota- 
tion have been performed by many authors, see e.g. 

EElESEol^^ 

Since the experiments by Rossby in 1969 |43], it is 
known that rotation can also enhance the heat transport. 
Rossby found that, when water is used as the convective 
fluid, the heat transport first increases when the rota- 
tion rate is increased. He measured an increase of about 
10%. This increase is counterintuitive as the stability 
analysis of Chandrasekhar [27 ] has shown that rotation 
delays the onset to convection and from this one would 
expect that the heat transport decreases. The mecha- 
nism responsible for this heat transport enhancement is 
Ekman pumping 03 H B3 |1B ED, i.e. due to 
the rotation, rising or falling plumes of hot or cold fluid 
are stretched into vertically-aligned-vortices that suck 
fluid out of the thermal boundary layers adjacent to the 
bottom and top plates. This process contributes to the 
vertical heat flux. For stronger rotation Rossby found, 
as expected, a strong heat transport reduction, due to 
the suppression of the vertical velocity fluctuations by 
the rotation. This means that a typical measurement of 
the heat transport enhancement with respect to the non- 
rotating case as function of the rotation rate looks like 
the one shown in figure [2] After Rossby many experi- 
ments, e.g. [f33] [26] [45] [23] [46] [20] [2T1 ITO] |T2l |T7] 
[47l|48l|4E[l5l[l6l[5a[5ll, have confirmed this general 
picture. A new unexplained feature observed in high 
precision heat transport measurements is the small de- 
crease in the heat transport just before the strong heat 
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Figure 2: The scaled heat transfer Nu(l/Ro)/Nu(0) as function of 
1/Ro on a logarithmic scale, a) Experimental and numerical data for 
Ra = 2.73 x 10 8 and Pr = 6.26 in a Y = 1 sample are indicated by 
red dots and open squares, respectively. Data taken from Stevens et 
al fl2l . 
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Figure 3: The Nusselt number Nu(l/Ro), normalized by Nu(0) with- 
out rotation, as a function of the inverse Rossby number 1 /Ro. Data 
are from experiments of Zhong and Ahlers 1 17 1. Main figure: r = 1, 
Pr = 4.38, Ra = 2.25 x 10 9 (solid circles), 8.97 x 10 9 (open circles), 
and 1.79 x 10 10 (solid squares). Inset: zoom in of figure [5] Figure 
taken from Weiss et al. fl3l . 



transport enhancement sets in, see figure [3] The experi- 
ments reveal that this effect becomes stronger for higher 
Ra. Because only for Ra > 10 10 the effect is more than 
1% it is currently not possible to study this in simula- 
tions, as at these high Ra numbers the spatial and time 
resolution would be insufficient to capture such a small 
effect. 

In early numerical simulations of rotating RB con- 
vection a horizontally unbounded domain was simu- 
lated. These simulations are relevant to understand 
large convective systems that are influenced by rotation, 



3 





• Stevens et al. sim (201 1 - 2009) 

• Stevens et al. exp (201 1 ) 

• Weiss etal. (2011 &2010) 

• Zhongetal. (2010 & 2009) 

• Niemela et al. (2010) 

• Tilgneretal. (2010 & 2009) 

• King et al. exp (2009) 



♦ King et al. sim (2009) 

♦ Liu and Ecke (2009 & 1997) 

♦ Kunnen et al. (2008) 
Oresta et al. (2007) 

♦ Kunnen et al. (2006) 

♦ Julien etal. (1996) 

♦ Rossby (1969) 



Figure 1: Phase diagram in Ra-Pr-Ro space for rotating RB convection; an updated version from Stevens et al. fill . The data points indicate 
where Nu has been measured or numerically calculated. The data are obtained in a cylindrical cell with aspect ratio r = 1 with no slip boundary 
conditions, unless mentioned otherwise. The data from DNSs and experiments are indicated by diamonds and dots, respectively. The data sets are 
from: Stevens etal (2011 - 2009) 1 10 12]|]2][l4][TT]|I] (T = 0.5, r = 1.0, T = 4/3, T = 2.0); Weiss & Ahlers (2011 & 2010) [[3[l5][l6) (r = 0.5) 
; Zhong & Ahlers (2010 & 2009) (K)l[l7]; Niemela et al (2010) (H (r = 0.5); Schmitz & Tilgner (2009 & 2010) 02 (free slip and no-slip 
boundary conditions and horizontally periodic); King et al (2009) (horizontally periodic and different aspect ratios) f20l ; Liu & Ecke (2009 & 
1997) E][22l (square with r = 0.78); Kunnen etal (2008) |23 1; Oresta etal (2007) f24l (T = 0.5); Kunnen etal (2006) |25 1 (r = 2, horizontally 
periodic); Julien et al (1996) | 26 | (r = 2, horizontally periodic); Rossby (1969) (varying aspect ratio). Panel a shows a three-dimensional view 
on the phase space (see also the movie in the supplementary material), b) projection on the Ra-Pr phase space, c) projection on the Ra-Ro phase 
space, d) projection on the Pr-Ro phase space. 
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such as the atmosphere and solar convection. In 1991 
Raasch & Etling [52] used a large-eddy simulation to 
simulate the atmospheric boundary layer. Julien and 
coworkers |[26j HU [53l EU, Husain et al (53, Kun- 
nen et al EU EH and King et al (20l [S3 used DNSs 
at several Ra and Ro to study the heat transport and 
the resulting flow patterns under the influence of rota- 
tion. Sprague et al l57l numerically solved an asymp- 
totically reduced set of equations valid in the limit of 
strong rotation. Finally, Schmitz & Tilgner (191 [58l 
performed horizontally periodic simulations with no- 
slip and stress-free boundary conditions at the horizon- 
tal plates. They show that heat transport enhancement 
under the influence of rotation is only found when a 
no- slip boundary condition at the horizontal plates is 
used. In recent years the focus of numerical simulations 
has shifted from horizontally periodic simulations to ro- 
tating RB convection in a cylindrical sample in order 
to make a direct comparison with experiments possible 
Ga[23|47l|48l|4^ As all 

information on the flow field is available in these simu- 
lations these studies focus on determining the influence 
of rotation on the heat transport and the corresponding 
changes in the flow structure. 

A direct comparison between experiments and simu- 
lations was used by Zhong et al ifTUl and Stevens et al 
lfT4l to study the influence of Ra and Pr on the effect 
of Ekman pumping. A strong heat transport enhance- 
ment with respect to the non-rotating case was observed 
for Pr « 6 IfTUl . The maximum enhancement decreases 
with increasing Ra and decreasing Pr, see figure [4] and 
[5] Later Stevens et al lfT4l found that at a fixed Ro 
number the effect of Ekman pumping, and thus the ob- 
served heat transport enhancement with respect to the 
non-rotating case, is highest at an intermediate Prandtl 
number, see figure [5] At lower Pr the effect of Ek- 
man pumping is reduced as more hot fluid that enters 
the vortices at the base spreads out in the middle of the 
sample due to the large thermal diffusivity of the fluid. 
At higher Pr the thermal boundary layer becomes thin- 
ner with respect to the kinetic boundary layer, where 
the base of the vortices is formed, and hence the tem- 
perature of the fluid that enters the vortices becomes 
lower. In this framework the decrease of the heat trans- 
port enhancement at higher Ra can be explained by the 
increase of the turbulent viscosity with increasing Ra. 
These observations also explain why no heat transport 
enhancement due to rotation is observed in the very high 
Ra number experiments of Niemela et al (THJ. In re- 
cent simulations of Horn et al 1 60] the effect of non- 
Oberbeck-Boussinesq (NOB) corrections was studied in 
rotating RB convection and they showed that for water 
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Figure 4: The ratio Nu(l/Ro)/Nu(0) as function of l/Ro for different 
Ra in a T = 1 sample. The experimental results for Ra = 2.99 x 10 8 
(Stevens et al QD), Ra = 5.63 x 10 8 (Zhong & Ahlers |17|), and 
Ra = 1.13 x 10 9 (Zhong & Ahlers 1 17]) are indicated in black, red, 
and blue solid circles, respectively. The DNS results from Stevens et 
al for Ra = 2.91 x 10 8 , and Ra = 5.80 x 10 8 are indicated by 
black and red open squares, respectively. All presented data in this 
figure are for Pr = 4.38. Figure taken from Stevens et al (8). 
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Figure 5: The heat transfer as function of Pr on a logarithmic scale 
for Ra = 1 x 10 8 . Black, red, blue, and green indicate the results for 
Ro = oo, Ro = 1.0, Ro = 0.3, and Ro = 0.1, respectively. The data 
obtained on the 385x 193x385 (open squares) and the 257x 129x257 
grid (solid circles) are in very good agreement. Figure taken from 
Stevens et al fTH 

NOB corrections leads to a Nusselt number that is a few 
percent higher than for the Oberbeck-Boussinesq case. 

4. Different turbulent states 

When the heat transport enhancement as function 
of the rotation rate is considered, a division in three 
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Figure 6: Visualization, based on a DNS at Ra = 2.91 x 10 8 with Pr = 4.38 in a T = 2 sample, of the temperature field in the horizontal 
mid-plane of a cylindrical convection cell. The red and blue areas indicate warm and cold fluid, respectively. The left image is at a slow rotation 
rate (l/Ro = 0.046) below the transition where the warm upflow (red) and cold downflow (blue) reveal the existence of a single convection roll 
superimposed upon turbulent fluctuations. The right image is at a somewhat larger rotation rate (l/Ro = 0.6) and above the transition, where 
vertically-aligned vortices cause disorder on a smaller length scale. 

regimes is possible l6T1l47ll49l . Here we will call these 
regimes: regime I (weak rotation), regime II (moder- 
ate rotation), and regime III (strong rotation), see figure 
[2] It is well known that without rotation a large scale 
circulation (LSC) is the dominant flow structure in RB 
convection, see e.g. Ref. [ 1 ] and figure|6^. This has mo- 
tivated Hart, Kittelman & Ohlsen 162IL Kunnen et al. 
l23ll . and Ahlers and coworkers (63l [171 [16) to study 
the influence of weak rotation on the LSC. In these in- 
vestigations it was found that the LSC has a weak pre- 
cession under the influence of rotation. In addition, it 
was shown by Zhong & Ahlers 1 17 ] that, although the 
Nusselt number is nearly unchanged in regime I (see 
figure [2]), there are various properties of the LSC that 
are changing significantly when the rotation rate is in- 
creased in regime I. Here we mention the increase in 
the temperature amplitude of the LSC, the LSC preces- 
sion (also observed by Hart et al. l62l and Kunnen et 
al. (23)), the decrease of the temperature gradient along 
the sidewall, and the increased frequency of cessations. 
To our knowledge these observations are still not fully 
understood. 

It was shown by Stevens et al. (12 ] that the heat trans- 
port enhancement at the start of regime II sets in as a 
sharp transition, see inset figure [2| due to the transition 



from a flow state dominated by a LSC at no or weak ro- 
tation to a regime dominated by vertically-aligned vor- 
tices after the transition. They demonstrate that in ex- 
periments the transition between the two states is in- 
dicated by changes in the time averaged LSC ampli- 
tudes, i.e. the average amplitude of the cosine fit to the 
azimuthal temperature profile measured with probes in 
the sidewall, and the temperature gradient at the side- 
wall. Later experiments and simulations by Weiss & 
Ahlers lfT6l and Kunnen et al. l49l confirmed that also 
other statistical quantities of the sidewall measurements 
change at the onset of heat transport enhancement. In 
addition, the simulations by Stevens et al. [12] and 
Weiss et al. fT3l show a strong increase in the number of 
vortices at the thermal boundary layer height, when the 
heat transport enhancement sets in. It was also revealed 
that at this point the character of the kinetic boundary 
layer changes from a Prandtl-Blasius boundary layer at 
no or weak rotation to an Ekman boundary layer af- 
ter the onset, which is revealed by the Ro 112 scaling of 
the kinetic boundary layer thickness after the onset, see 
figure [7] Furthermore, a strong increase in the vertical 
velocity fluctuations at the edge of the thermal bound- 
ary layer, due to the effect of Ekman pumping, was ob- 
served, while the volume averaged value decreased due 
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to the destruction of the LSC. 

Here we also mention that there have been analytic 
modeling efforts by Julien et al. (64), Portegies et al. 
lf65l and Grooms et al. [ 66 ] to understand the heat trans- 
port in regime III. In these models only the heat transfer 
in the vertically-aligned vortices is considered as most 
heat transport in this regime takes place inside these vor- 
tices. The trends shown by these models are in good 
agreement with simulation results. Furthermore, it was 
shown by Kunnen et al. [ 49 ] that under the influence of 
rotation a secondary flow is created in regime II and III. 
This flow is driven by the Ekman boundary layers near 
the plates and generates a recirculation in the Stewart- 
son boundary layer on the sidewall with upward (down- 
ward) transport of hot (cold) fluid close to the sidewall 
in the bottom (top) part of the cell. Hence this secondary 
flow leads to the generation of a vertical temperature 
gradient along the sidewall, which is measured in exper- 
iments |49l EOl US [13 and simulations [49]. We note 
that this secondary circulation is significantly different 
from that suggested by Hart et al. (67 ], based on a linear 
mean temperature gradient, or that of Homsy & Hudson 
[68]. Under the influence of strong rotation, i.e. for 
small values of Ro, a destabilizing temperature gradi- 
ent is also formed in the bulk. This temperature gra- 
dient is caused by the merger of the vertically-aligned 
plumes [see, e.g.,[69l[33j|26l|70l[53l[54^ i.e. the en- 
hanced horizontal mixing of the temperature anomaly of 
the plumes results in a mean temperature gradient in the 
bulk. Recently, this effect was measured for different 
1 /Ro and Ra by Liu & Ecke [ 71 ] with local temperature 
measurements in the bulk, which show good agreement 
with earlier simulations [10 ] and experiments by Hart et 
al. (671. 



5. Sidewall temperature measurements 

In recent high precision heat transport measurements 
of rotating RB convection the samples are equipped 
with 24 thermistors embedded in the sidewall (lOl [T71 
[T3l[T6lll9]|. These thermistors are divided over 3 rings 
of 8 thermistors that are placed at 0.25z/L, 0.50z/L, and 
0.75z/L. In non-rotating convection this arrangement of 
thermistors is used to determine the orientation of the 
LSC as the thermistors can detect the location of the up- 
flow (downflow) by registering a relatively high (low) 
temperature. This method has been very successful in 
studying the statistics of the LSC in non-rotating RB 
convection, see the review of Ahlers et al. Q. 

To determine the strength and orientation of the LSC 
a cosine function is fitted to the azimuthal temperature 
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Figure 7: The thickness of the kinetic top (diamonds) and bottom 
(circles) boundary layer for Ra = 2.73 x 10 8 and Pr = 6.26 in a T = 1 
sample. The black line indicates the Ekman BL scaling of (l/Ro)' 1 ^ 2 . 
Figure taken from Stevens et al. fTTl . 



profile obtained from the sidewall measurements (T) as 

6i = k + 5 k cos(0/ - cf) k ). (11) 

Here 0; = 2in/n p , where n p is the number of probes, 
indicates the azimuthal position of the probes, while k 
gives the mean temperature, S k the temperature ampli- 
tude of the LSC, and (p k the orientation. The index k 
indicates whether the data are considered at z/L - 0.25 
(b), z/L = 0.50 (m) or z/L = 0.75 (t). In simulations 
similar measurements are obtained by Stevens et al. (59l 
and Kunnen etal. (23l[49) from numerical probes placed 
close to the sidewall. 

In rotating RB convection it is important to know 
when the LSC breaks down. Zhong et al. [10] showed 
that at the onset of heat transport enhancement the tem- 
perature amplitude of the LSC decreases, while an in- 
crease in the vertical temperature gradient at the side- 
wall is observed. Kunnen et al. (23ll used simulation 
data to show that the energy in the first Fourier mode 
of the azimuthal temperature profile at the sidewall de- 
creases strongly when the heat transport enhancement 
sets in. This idea was quantified by Stevens et al. [72 ] 
who proposed to determine the relative LSC strength 
based on the energy in the different Fourier modes of 
the measured or computed azimuthal temperature pro- 
file at or nearby the sidewall, as 



S k = max 



t dt 



1 

N 



I 1- 



,0 



(12) 



The index k has the same meaning as indicated below 
eq. |lTJ). In eq. |l2| E\dt indicates the sum of the 



energy in the first Fourier mode over the time interval 
[tb,t e ] of me simulation or experiment, and f f e E tot dt 
the sum of the total energy in all Fourier mocfes over 
the same time interval. N indicated the total number of 
Fourier modes that can be determined from the number 
of azimuthal probes that is available. From the defini- 
tion of the relative LSC strength it follows that always 
< Sk < 1. Concerning the limiting values: Sk = 1 
means the presence of a pure azimuthal cosine profile 
and Sk = indicates that the magnitude of the cosine 
mode is equal to (or weaker than) the value expected 
from a random noise signal. We consider a value for S k 
of about 0.5 or higher as an indicated that a cosine fit 
on average is a reasonable approximation of the data set 
and Sk well below 0.5 that most energy is in the higher 
Fourier modes. It was shown by Kunnen et al. |49], 
Stevens et al. |59], and Weiss et al. lfT6l that in a T = 1 
sample the relative LSC strength is close to 1 before the 
onset of heat transport enhancement and strongly de- 
creases to values close to zero at the moment that heat 
transport enhancement sets in, thus confirming that the 
LSC breaks down at this moment, see figure [8^. 

6. Determination of vortex statistics 

In regime II and III, see figure [2j the flow is dom- 
inated by vertically-aligned vortices. The experiments 
of Boubnov & Golitsyn (6T|, Zhong, Ecke & Steinberg 
l33lL and Sakai 1731 first showed with flow visualization 
experiments that there is a typical ordering of vertically- 
aligned vortices under the influence of rotation, see fig- 
ure [6] for a visualization based on numerical simula- 
tions. A three-dimensional visualization of the vortices 
is made by Kunnen et al. l48l and Stevens et al. (59), 
also based on numerical simulations. First experimental 
measurements of flow and temperature in rotating RB 
samples have already been reported some 20 years ago. 
It started with local velocity (and temperature) measure- 
ments in the flow have been reported by Boubnov & 
Golitsyn ifSTIl and Fernando et al. [ 74 1. In both ex- 
periments the flow was analyzed using particle- streak 
photography. Later Vorobieff & Ecke El |45) used 
digital cameras to perform particle image velocimetry 
(PIV) measurements to study the flow with higher spa- 
tial and temporal resolutions. In addition, they added 
thermochromic liquid crystal (TLC) microcapsules to 
visualize the vortices in the temperature field. In later 
experiments by Kunnen et al. [23 , 46 , 47] flow measure- 
ments were extended by the application of stereoscopic 
PIV (SPIV): In this technique two cameras are used to 
obtain a stereoscopic view, which allows for a recon- 
struction of the out-of-plane velocity. This made new 



experimental data available, for example on the corre- 
lation between the vertical velocity and vorticity. Re- 
cently, also King et al. 1 76 1 used flow visualization ex- 
periments to confirm the formation of vortices in the ro- 
tating regime. 

Although the presence of vortices in many cases is 
clear from snapshots of the flow, a well-defined vor- 
tex detection criteria are necessary to study the vortex 
statistics. Common methods are based on the velocity 
gradient tensor Vu = diUj(iJ <E 1,2,3) 177117811 . This 
tensor can be split into a symmetric and antisymmetric 
part 

Vu = i[Vu + (Vu) r ] + i[Vu - (Vu) r ] = S + D. (13) 
The Q^d criterion 17711 defines a region as vortex when 

Q3D= i(||n|| 2 -||^|| 2 )>0, (14) 

where ||A|| = ^Tr(AA T ) represents the Euclidean norm 
of the tensor A. The two-dimensional equivalent is 
based on the two-dimensional horizontal velocity field 
perpendicular to z, i.e. on Vu^d = diUj(i,j e 1,2), 
which defines a region as vortex when 1751 l45l 

Q 2D ee 4Det(Vu| 2D ) - [7>(Vu| 2D )] 2 > 0. (15) 

The Q2D criterion is used by Stevens et al. (T21 El and 
Weiss et al. lfT3l to determine the vortex distribution 
close to the horizontal plate. This analysis has con- 
firmed that vortices are indeed formed close to the hor- 
izontal plates. In addition, the analysis revealed that no 
vortices are formed close to the sidewalls while their 
distribution is approximately uniform in the rest of the 
cell. It was shown by Kunnen et al. l48l that the Qid cri- 
terion is not suitable to recover the smaller scale vortices 
in the bulk, where the Q^d criterion is more suitable. 
In addition, they present statistical data, based on sim- 
ulations and SPIV measurements, on the number and 
size of vortices at different locations in the flow. Later 
Stevens et al. l59l compared the regions that are iden- 
tified as vortex using the Q^d criterion with the regions 
of warm (cold) fluid that are found with temperature iso- 
surfaces. They show that the vortex regions found with 
the Qid criterion agree well with the regions shown by 
warm (cold) temperature isosurfaces, see figure [9] The 
main difference is that smaller vortices in the bulk are 
not shown by the temperature isosurfaces, because their 
base is not close to the bottom (top) plate where warm 
(cold) fluid enters the vortices. 
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Figure 8: The relative LSC strength at 0.50z/L, i.e. S m , as function of 1/Ro based on experimental data (open squares) and simulations (filled 
circles). Panel (a) shows experimental data from Zhong & Ahlers fT71 at Ra = 2.25 x 10 9 and simulation results from Stevens et al. l59l at 
Ra = 2.91 x 10 8 in a T = 1 sample. Panel (b) shows experimental data from Weiss & Ahlers [T6) and simulation results from Stevens at al. (59] 
for Ra = 4.52 x 10 9 in a Y = 1/2 sample. All presented data are for Pr = 4.38. The vertical dashed line indicates the position where heat transport 
enhancement starts to set in, i.e. at l/Ro c w 0.86 (l/Ro c « 0.40) for r = 1/2 (r = 1) according to Weiss et al. fT3l . Figure taken from Stevens et 
al (59) 





Figure 9: Flow visualization for Ra = 4.52 x 10 9 , Pr = 4.38, and 
l/Ro = 3.33 in a T = 1/2 sample. The left image shows a three- 
dimensional temperature isosurfaces where the red and blue region 
indicate warm and cold fluid, respectively. The right image shows a 
visualization of the vortices, based on the Q^d criterion. Note that the 
red (blue) regions indicated by the temperature isosurface correspond 
to vortex regions. Figure taken from Stevens et al. (59l . 



wall by a strong decrease of the relative LSC strength, 
see figure [8^. However, this is not the case for all as- 
pect ratio samples. Namely, it was revealed by Weiss & 
Ahlers lfT6l that in a T = 1/2 sample no strong decrease 
in the relative LSC strength is observed at the moment 
that the heat transport enhancement sets in, see figure 
[3b. These authors discuss that the relative LSC strength 



according eq. ( 12) to does not allow one to distinguish 
between a single convection roll and a two- vortex- state, 
in which one vortex extends vertically from the bottom 
into the sample interior and brings up warm fluid, while 
another vortex transports cold fluid downwards. Be- 
cause a two- vortex- state results in a periodic azimuthal 
temperature variation close to the sidewall which can- 
not be distinguished from the temperature signature of 
a convection roll with up-flow and down-flow near the 
side wall. The consequence is that various other crite- 
ria, which are based on identifying a cosine variation of 
the temperature close to the sidewall, also do not show 
that a transition takes place. Stevens et al f59\ con- 
sidered this r = 1 /2 case in numerical simulations and 
they showed with flow visualization, vortex detection, 
and the analysis of several statistical quantities that in a 
r = 1 /2 sample indeed a flow state that on average can 
be considered as a two-vortex state is formed. 



7. Influence of the aspect ratio 

In a r = 1 sample the onset of heat transport enhance- 
ment is visible in temperature measurements at the side- 



The aspect ratio of the sample does not only deter- 
mine the number of vortices that is found, but also the 
rotation rate at which heat transport enhancement sets 
in, see figure [10] Weiss et al [Q21 Q21 showed with 
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Figure 10: The ratio Nu(l/Ro)/Nu(0) as function of l/Ro for /to « 
3 x 10 8 and different T. The experimental results for T = 1, T = 4/3, 
and r = 2 are indicated in red, dark green, and blue solid circles, 
respectively. DNS results for T = 0.5, r = 1, and r = 2 are indicated 
by black, red and blue open squares, respectively. All presented data 
are for Pr = 4.38. Figure taken from Stevens et al. (8). 



a Ginzburg-Landau-like model that the rotation rate at 
which the onset of heat transport enhancement sets in 
(l/Ro c ) increases with decreasing aspect ratio due to fi- 
nite size effects. Based on experimental and numerical 
data they find that heat transport enhancement sets at at 
critical Rossby number Ro c 



1 

RoZ 



at b 

r 1 + f 



(16) 



where a = 0.381 and b = 0.061. In addition, predictions 
of the theory about the horizontal vortex distribution in 
the region close to the sidewall (131 [15| are closely re- 
produced by numerical measurements (T3HH1 - Although 
the aspect ratio is important for the heat transport at rel- 
atively weak rotation it is shown by Stevens et al. 
that the aspect ratio dependence of the heat transport 
disappears for sufficiently strong rotation rates. 

8. Conclusions 

We have summarized and discussed our work on ro- 
tating Rayleigh-Benard (RB) convection over the last 
years and how it is related to other work on rotating RB. 
We have seen that a combination of experimental, nu- 
merical, and theoretical work has greatly increased our 
understanding of this problem. As is shown in the rotat- 
ing RB phase diagrams, see figure [T] some parts of the 
phase diagram are still relatively unexplored. We es- 
pecially note that up to now the rotating high Ra num- 
ber regime has only been achieved in the experiments 



of Niemela et al. lfT8l . Because many natural phenom- 
ena like the convection in the atmosphere 1 4 ] and oceans 
0, including the global thermohaline circulation (6j, 
are influenced by rotation. In all these cases the Ra 
number is very high, and it is thus very important to 
understand this regime better. It is especially important 
to know how rotation influences the different turbulent 
states that are observed in the high Ra number regime 
l79l . In addition, it would be very nice to have more 
detailed data for high and low Pr numbers over a much 
larger Ra number range, as nowadays most datasets are 
clustered around Pr = 4.3 and Pr = 0.7. Furthermore, 
there are still several issues that need a deeper under- 
standing. Here we mention the relatively strong changes 
in the LSC characteristics that are observed in regime I 
and the small decrease in the Nusselt number observed 
in high Ra number experiments just before the onset of 
heat transport enhancement. It would also be very inter- 
esting to get more detailed information about the statis- 
tical behavior of the vortical structures in regime II and 
III and to see whether it is possible to describe all heat 
transfer properties in rotating RB convection by a uni- 
fying model like in the Grossmann-Lohse theory [ 1 ] for 
non-rotating RB convection. 
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